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I . INTRODUCTION 


The  atmospheric  detonation  of  a nuclear  weapon  produces  a high 
temperature  air  plasma  which  is  commonly  referred  to  as  a "fireball". 

This  fireball  begins  to  emit  thermal  radiation,  mostly  in  the  visible 
and  infrared  regions  of  the  spectrum,  into  the  surrounding  atmosphere 
immediately  after  it  is  formed.  Although  the  fireball  has  a finite 
diameter,  it  is  considered  to  be  a point  source  for  calculational  purposes, 
and  the  thermal  radiation  environment  produced  is  expressed  in  terms  of 
a time  dependent  thermal  irradiance.  The  general  characteristics  of  this 
irradiance  at  any  point  outside  the  fireball  are  shown  by  the  curve  in 
Figure  1.  Jhis  curve  }s  generally  approximated  by  2 


H(t)  - H, 


t j pvxu v 

curve  i 


0.206  e 


1 ♦ 1.69 


W 


(1) 


where  t is  a function  of  the  weapon  yield  and  Hq  is  a function  of 
weapon  yield,  distance  from  point  of  weapon  detonation,  and^transmittance 
of  the  atmosphere.  Table  1 contains  values  of  tQ,  Hq  andQ/  H(t)  dt  for 
several  weapon  yields  and  distances.  In  theory  the  thermal  environment 
is  of  infinite  duration. 


TABLE  I.  Nuclear  Pulse  Parameters 


Yield 

(Kt) 

Distance 

(km) 

(s? 

Ho  , 
(MWm  ) 

PO 

J Hd*2 

°(MJm  *) 

1 

0.32 

0.04 

11.32 

1.25 

10 

0.70 

0.11 

11.26 

3.26 

100 

1.51 

0.31 

8.96 

7.31 

1000 

3.26 

0.87 

4.65 

14.08 

However,  for  practical  applications  the  time  duration  of  the  environment 

lOt 


is  taken  to  be  10tQ  since 


/ 


H(t)dt> 


0.85  J 


H(t)dt  and  H(t)<0.03  H 


for  t > lOt  . 

o 


The  effects  of  the  thermal  environment  on  exposed  materials  are  due 
to  the  absorption  of  all  or  part  of  the  radiant  energy  incident  on  the 
exposed  surfaces  of  the  materials.  The  absorption  of  this  radiant  energy 
by  the  lateral  surface  of  a cylinder  will  result  in  a transient  temperature 
rise,  a knowledge  of  which  is  essential  for  predicting  the  response  of 
the  cylinder  for  those  effects  which  are  a function  of  temperature  or 
temperature  change.  Because  of  the  transient  heating,  a numerical 
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Figure  1.  Nuclear  Thermal  Pulse 


calculation  of  the  temperature  using  finite  difference  or  finite  element 
methods  is  generally  employed.  Although  these  methods  are  of  great 
value  for  solving  a specific  problem,  they  are  not  convenient  as  analytical 
methods  for  nondimensional  representation  or  parametric  analysis  of  the 
temperature  field  in  the  cylinder. 

This  report  describes  a deviation  of  the  transient  temperature  field 
in  an  isotropic,  homogeneous,  finite  length,  solid  cylinder  of  radius  r 
whose  lateral  surface  is  subjected  to  heating  by  a nuclear  thermal 
radiation  environment.  Ojalvo's  modified  separation-of-variables  method1 
is  used  to  solve  the  transient  heat  conduction  equation  under  the 
following  assumptions: 

1.  the  thermal  properties  of  the  cylinder  are  independent  of 
temperature, 

2.  the  thermal  radiation  is  absorbed  at  the  surface, 

3.  convection  and  radiation  heat  losses  by  the  cylinder  can  be 
neglected, 

4.  the  transient  temperature  field  in  the  cylinder  is  independent 
of  the  axial  coordinate,  and 

5.  the  initial  temperature  field  in  the  cylinder  is  uniform. 


II.  TEMPERATURE  EQUATIONS 


The  thermal  irradianee  at  the  lateral  surface  of 
expressed  by 


H(6,t)  « 


) cos  0 

0 


O£0 

2s”  2 


the  cylinder  is 


(2) 


j i*0f(t)  cos  0 3n  so  £2n 

where  f(t)  is  the  time  dependent  portion  of  (1).  The  transient  temperature 
fiold  in  the  cylinder  is  governed  by  the  following  equations 


3T(r,6t) 

0t 


(3) 


whore 


7 - t 3?  Cr  3r  > * r2  JP  * 


(4) 


1 J.U.  Qjalvo,  "Conduction  with  Tim-Dependent  Heat  Sources  and  Boundary 
Conditions,"  International  Jouintal  of  Heat  and  !4ao8  Transfer,  Vol  5, 
1962,  pp.  1106-1109. 
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3T 

v9r 


JL,f  cos  8 o<e<|- 


K 

0 


(5) 


Hpf  cos  8 3n  < 8 < 211 

K 2 


at  r = ro> 


1 9T 

r 38  “ °* 

at  8 = 0,  n (because  of  symmetry),  and 

T(r,8,0)  = 0*. 


(6) 


(7) 


The  modified  separation-of-variables  method  assumes  that  the 
solution  of  (3)  can  be  expressed  as 


T 1 E E *«ntt)*«»tr'W  * ToCr,0)«t). 


(8) 


m n 

By  substituting  (8)  into  (3)  one  has 


EE*  *2*  * * t *1 

Tmn  vmn  o < / Jm  ^mn  o I . (9) 


The  method  further  assumes  that 


V""* 

T ■ > a $ 
o / j am  on 


and  that  either 


VT  2 » 0 
o 


or 


vh  » y*  b <J> 

o / j anwon 


(10) 


(11) 


(12) 


So  generality  ia  loot  by  as  aiming  T * 0 since  the  temperature  of  the 
cylinder  ocn  be  obtained  by  adding  the  initial  constant  tenpeiviture 
to  T and  only  the  temperature  difference  is  required  for  the  sires* 
and  displacement  fields. 


The  substitution  of  (10)  and  (12)*  into  (9)  yields 


V~  4> 
cm 


<!> 


mn  (i3) 


If  (13)  is  equated  termwise  and  divided  by  and  if  - y's  are 

chosen  for  the  separation  constants,  one  obtains  the  following 
differential  equations: 


and 


V2  db  ♦ X2  <i>  =0 

*mn  mn  *mn 


i +X2  — * = -a  f + — b f 

Tmn  m cp  u mn  c.p  mn 


(14) 

(15) 


The  solution  of  (14)  and  (15)  along  with  Tq  when  substituted  into  (8) 
provides  the  transient  temperature  field  in  the  cylinder. 


The  boundary  and  initial  conditions  for  <Jt>  , T , 
obtained  by  substituting  (8)  into  (5),  (61,  and  (7). 
yieu,s 


EE* 


3$ 


on 


mn  3 r 


3Tf 
3r 


at  r * rQ, 


III  cos  8 2T„  e 

“O  - r~~0  I 

«■  3r 

m • 


V > V ’A 

LuLa 


* . Jla.  f 

an  r38  r38 


n 

T 

3ft 


7 

<0<~ 

w 

<0  <28 


and  ♦ are 
ran 

This  substitution 


(16) 


(17) 


at 8 ■ 0,11,  and 


EE*-.  4u  • -v 


(18) 


rhe  use  of  (22)  is  mctfw  general  than  the  use  of  (11)  since  ill)  is  a 
special  case  Cf  (22)  in  which  all  the  oiv  zwo. 
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at  t = 0.  From  (16),  (17),  (18)  and  Figure  1 the  following  boundary  and 

initial  conditions  for*  , T , and  ♦ are  deduced: 

mn  o mn 


at  r = r , 
o 


at  = Q,  3 , n 


3T 

3r 


3^mn  a 


Tr 


34- 


mn 


r80 


at  r = r , 
o’ 


at  3 o o,  "• , and 
at  t * 0. 


Hp  cos  8 

K 

0 

H 

-o  cos  0 

K 


rSW°  * 0 


♦mn  C 0 

on 


o < e<n 

2 

— $0<y- 

~~<o<2n 


(19) 

(20) 

(21) 

(22) 

(23) 


liquation  (14)  is  a two  dimensional  Helmholtz’s  differential 
equation,  the  general  solutions  of  which  are  of  the  form 

**>«,  ■ Jo  <V>  Uo»  * Bon0)  (24“) 

for  o *!  0 and 

$®i  “ {AmnJaUmnr)  * V(Xm,»r)  HCacos  ffl0  * Vin  ffl0)  (24b) 

for  tax). 

In  order  for  these  solutions  to  be  well  behaved  at  r «*  0,  theB  ’s  must 
be  equal  to  tero.  Also  from  (20)  it  is  readily  seen  that  the  Bon's  and 
the  Dg's  must  also  equal  zero  and  that  a » 0,  1,  2,  3,  , Equation 

(24)  can  now  be  rewritten  as  * 


*)  COS  a8 

an  on  n tan 


(2S) 


This  equation  ineupea  a kaaogeneoua  bowtdapy  vendition  for  (14) . 
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The  substitution  of  (25)  into  (19)  gives 


or 


J (X  r ) o 0 
m ran  o 


rx  r 'i  Ja>l^mnro^  _ 

on  o J (X  r l 
m mn  o 


(26a) 

(26b) 


From  (26a)  one  sees  that  the  * r 's  ar©  the  n-th  positive  roots  of 
(26a) . 

Equation  (11)  cannot  be  used  in  determining  T . The  solution  of 
this  two  dimensional  Laplace’s  differential  equation  which  is  well 
behaved  at  r = 0 and  satisfies  (22)  is 

oo 

T = C ♦ V C rm  cos  n 0 (27) 

o o f j m 

o=l 

the  substitution  of  (27)  into  (21)  yields 

H O<0^4~ 

— cos  0 2 

1C 

o !<o£-jr  (2S) 

- cos  0 ^-<9  < 211 

K 2 


E 

o»l 


aC  r 
a o 


m-1 


COS  00s 


By  expressing  the  right  side  of 
cosine  series: 


«o 


i 

n 


t '1  9 
♦ «-*  - <■— 


(23)  in  terms  of  the  following  Fourier 


n 


tni 

/ j. 


E w-n 

n«2,4,6,, 


cos  n 0 


(29) 


one  can 

50 


ffl“l 


rewrite  (28)  as 

aC  r cos  a 9 
b o 


n 


iio 


l 

rf 


L 


. COS  9 
* 


even 


it  is  readily  seen  from  (30)  that  (27)  cannot  satisfy  the  boundary 
conditions  at  f « f because  of  the  first  term  of  (27).  This  term 
represents  the  uniform  heating  of  the  entire  lateral  surface  of  the 
cylinder,  and  the  form  of  T which  is  obtained  by  soiling  the  uniform 
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heating  case  is 


T = -ikL-r2 

o 2iltcr 


(31) 


Consequently,  it  is  assumed  that 


oo 


T « 
o 


H 


m 


2IItcr 


— 2 ♦ C + > C r cos  m 8 
r 0 ' 4 m 

) m»l 


(32) 


where  the  values  of  the  C^'s,  except  for  C , are  obtained  by  equating 
(30)  termwise.  From  (30)  one  sees  that 


and 


C,  . 

1 7< 

(-) 

c = - 2Hn  cosy  2 / 
m nm(m2-l)icron" 1 


for  o > 1 . 

In  order  to  solve 


(33) 

(34) 


• k*  2 * ic 

y « — \c  q>  o -a  f ♦ l>  f 

nm  cp  mn  mn  an  cp  an 


(15) 


the  coefficients,  a_„  and  b_,  . oust  bo  evaluated.  The  a *s  can  be 

ffij)  rail  (till 

evaluated  from  (101  and  (32)  and  the  fc  's  from  (12)  and  (32).  From 
(10)  and  (32) 


tan 


oo  oo 

Jo-r2  . C .Vc  .vy.  j (X  r)  e 

2a«ro  o JLmJ  a oo  / -< / ^ tan  o'  tan 

»«0  n*l 

Equating  (35)  termwise  one  obtains 


cos  m 8 
(35) 


if  _2 


2Hx-r_ 


and 


♦ C * a ♦ 
o oo 


C r 
u 


n»l 


7 a „ d (*  ”) 
/ ri  on  o'  on  * 


a*E 


n»l 


a J (X  r). 
an  n an 


(36) 


(37) 
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Using  the  orthogonality  properties  of  the  Bessel  functions,  the  a 's 
can  be  evaluated  from  m 


Jki_ 

2Mv 


So  So 

J r3jo(Xonridr  + Co  JTj 


(X  r)dr 
o on  J 


(38) 


on 


A 


C\mr)dr 


and 


/ r”*1 


a = C 
mn  ra 


J (X  r) dr 
m mn  1 


J 


rJ  (X  r)dr 
m v mn  J 


The  integi'ation  of  (38)  yields 


a = + c 

oo  n<  o 


(39) 


(40p) 


for  X =0  and 
oo 


iijn£a. 


on  IIk(X  r )2J  (X  r) 
on  o'  o on  ' 


(40b) 


for  n*l.  For  che  remaining  a^'s  the  integration  of  (39)  yiolds 


m+1 


a 3 
mn 


20  X r * J . (X  r ) 
m ran  o m+1  mn  o' 


f(X  v )2-n2ld2  (X  x- 
[ mn  o'  J m v mn  c 


) 


By  substituting  (26b) , (33)  and  (34)  into  (41)  one  obtains 

,2 


a . 3 
ml 


for  m>l  and 


a » 
ran 


* ^Nnlo } 2l]  J l(xmlro)~ 

2H0r2  cos  (mn/2) 

— r 2~tt m ““ 

Hk  nM  (X  rj -of  JM(Xmnr  ) 
v mn  q m mn  o 


for  ra>l  and  n>l.  From  (12)  and  (32) 
m»0  n-0 


(41) 


(42a) 


(42b) 


(43) 
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Using  the  orthogonality  properties  of  the  Bessel  functions  and  the  cosine 

functions,  the  b ’s  can  be  evaluated-from 

ran  r_  2ji 


ran 


2Hn  o 


/rVW>dr  / 


cos  ra6d0 


Iter^  rQ  211 

/rJ2(X  r )dr  / cos2m0d6 
m mn  o'  J 
o o 


(44) 


The  integration  of  (44)  yields 


2HC 


00 


Ikr 


(45a) 


for  m>0  and  n>  1. 


b = 0 
ran 


For  » 0 equation  (IS)  reduces  to 


(45b) 


dr  a ~a  f + ■ b f % 

*00  00  CP  OO  * 


(46) 


the  integration  of  which  is 


V , e - a f(t)  ♦ — « — b 
ol  oo  v * cp  oo 


j f(t)dt  ♦ 


00 


(47) 


From  the  initial  condition  equation  (23),  one  sees  that  B must  equal 
zero*  For  the  remaining  X ’s,  the  g9neral2solution  of  (iS)  can  be 
obtained  by  using  the  integrating  factor  X^  t 


The  use  of  this  factor  results  in 


-V 


« -a  e 
Ban  an 


tut  cp 


cp 


t 2 * 

A t 


t r ' X t - X — t . 

J omc°  f(t)dt  * B^e 


(46) 


One  also  sees  from  (23)  that  the  B^’s  oust  equal  zero.  The  integration 
of  the  right  hand  side  of  (48)  by  parts  yields 


♦ a*** 

an  an  w am  cp 


-X  t - 

s an  cp 


J x2  t 
* • » f(t)dt 


(49) 
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The  substitution  of  (25) , (32) , (40) , (42) , (45) , (47)  and  (49)  into  gives 


t.-Sb-I  f 

nfcro|  J 


Oa 


• X t 

on  cp 


f(t)dt  +i 


n»l  J (Xr) 
o on  o 


/ 


Xon  ^ 1 


xf(t)dt 


1 - “ r,\  t f 'xin  cp  1 

J.»„r>  * -!-£  -l  ■?"  


o'  on 


n»l 


B>1»V  - A Wo> 


A,_  ~~  t 


j e ln  CP  • 


f(t)dt 


OO  <»  r,  r I2 

V V'  lXm  o] 


cos(mn/2)  e 


-X  — — t 

mn  cp 


OWo)2"”2]  VWo> 


/ 


e>i  CP  f(tldt 


* Wr)  cos  ° 6 


(SO) 


whore  f(t)  is  the  time  dependent  portion  of  (!}. 


Ill*  NUMERICAL  RESULTS 

For  simplicity  and  generality,  the  derived  temperature  field 
equation  will  be  numerically  evaluated  in  terms  of  the  following 
dimensionless  quantities 


c 

r 


r 


6 - 


t 


* 

T 


0pcro 


(51) 


IS 


In  addition,  (1)  is  approximated  by  the  following  Fourier  series: 

1=1  0 


f(t)  = H 


Co 

2 


Ut 


- I Sin  It 


(52) 


If  (1)  were  used  in  evaluating  (50),  the  integrals  would  have  to  be 
computed  numerically.  Since  the  use  of  (52)  allowed  for  the  direct 
integration  of  these  integrals,  such  a representation  is  both  logical 
and  convenient,  regardless  of  the  number  of  terms  required  for  accuracy. 
All  of  the  calculations  were  carried  out  by  using  a IBIIVAC  1108  digital 
computer  and  the  number  of  terms  used  in  evaluating  the  double  series 
solution  were  such  as  to  insure  at  least  three  digit  convergence. 

Figures  2;  5,  and  4 are  plots  of  the  radial  temperature  distribution 
for  varior 1 values  of  0 , t and  8.  A.  comparison  of  the  temperatures 
calculated  usinc-  (50)  and  those  calculated  using  CINDA-3G2,  a finite 
differencing  he  t transfer  computer  program,  shows  a five  percent  or  less 
difference  in  the  calculated  temperatures.  Approximately  ninety  seconds 
of  machine  time  is  required^to  calculate  the  radial  temperatures  using 
(50)  for  psiied  values  of  $ and  t and  two  values  of  8. 


IV.  CONCLUSION 

The  analytical  expression  of  the  transient  temperature  field  in  an 
isotropic,  homogeneous,  finite  length,  solid  cylinder  whose  lateral 
surface  is  subjected  to  the  heating  by  a nuclear  thermal  radiation 
environment  has  bean  derived.  This  expression  provides  a convenient 
moans  for  the  nondimensionai  representation  and  parametric  analysis  of 
the  temperature  field  in  th“  cylinder.  In  addition,  this  temperature 
equation  can  be  used  in  tue  analysis  of  thoso  effects  dependent  on 
temperature  or  temperature  Jtange,  e.g.,  thermal  stress  in  a cylinder. 
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GLOSSARY  OF  TERMS 


mn 


ran 

c 

f.f(t) 

h 

r,9,z 

r 

r\t  ,T 

t 

t 

o 

m’  mn*  ra 
Bran*  Cw*  Dm 
Cl 
Cra 

H 

Jm[X)* 

V<*> 


T,T(r,9,t) 

Vo(r*9) 

8 

tc 

Non 


*m 


= Coefficients  defined  by  (38)  and  (39) 

= Coefficients  defined  by  (44) 

= Specific  heat 
= Time  dependent  portion  of  (1) 

= Half  length  of  cylinder 
= Cylindrical  coordinates 
= Radius  of  cylinder 

= Dimensionless  variables  defined  by  (51) 

= Time 

= Rise  time  of  nuclear  thermal  pulse 

* Unknown  coefficients  in  (24) 

= Coefficient  defined  by 
= Coefficient  defined  by 

= Time  dependent  irradiance  of  nuclear  thermal 
value 

= Maximum  irradiance  of  nuclear  thermal  pulse 
a Ordinary  Bessel  function  of  argument  x 

= ■»„(») 
dx 

» Temperature  in  cylinder 
= Part  of  solution  to 
a Dimensionless  variable  defined  by 

» Thermal  conductivity 
= Separation  constant  of  (13) 

* Density 

a Part  of  solution  to  (3) 

= Part  of  solution  to  (3) 
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